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ABSTRACT 

We study the build-up of the bimodal galaxy population using the NEWFIRM Medium-Band Survey, which 
provides excellent redshifts and well-sampled spectral energy distributions of k, 27, 000 galaxies with A' < 22.8 
at 0.4 < z < 2.2. We first show that star-forming galaxies and quiescent galaxies can be robustly separated with 
a two-color criterion over this entire redshift range. We then study the evolution of the number density and 
mass density of quiescent and star-forming galaxies, extending the results of the COMBO- 17, DEEP2, and 
other surveys to z = 2.2. The mass density of quiescent galaxies with M > 3 x lO'^M© increases by a factor of 
~ 1 from z ^ 2 to the present day, whereas the mass density in star-forming galaxies is flat or decreases over the 
same time period. Modest mass growth by a factor of ~ 2 of individual quiescent galaxies can explain roughly 
half of the strong density evolution at masses > IO^Mq, due to the steepness of the exponential tail of the 
mass function. The rest of the density evolution of massive, quiescent galaxies is likely due to transformation 
(e.g. quenching) of the massive star-forming population, a conclusion which is consistent with the density 
evolution we observe for the star-forming galaxies themselves, which is flat or decreasing with cosmic time. 
Modest mass growth does not explain the evolution of less massive quiescent galaxies 10'"^ Mq), which 
show a similarly steep increase in their number densities. The less massive quiescent galaxies are therefore 
continuously formed by transforming galaxies from the star-forming population. 
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L INTRODUCTION 

Large surveys, such as the Sloan Digital Sky Survey 
(SDSS), have begun to sample representative volumes of the 
nearby Universe. One of the more surprising results from 
these surveys is the existence of a bimodal galaxy population, 
manifested in correlations between a wide variety of galaxy 
properties both ob served (i.e., color vs. luminosity, color 
vs. morph ology; IStrateva et all 120011: IBlanton et all 120031; 
iBaldrv etal . 2004) and derived (i.e., stellar age and stellar 
mass; e.g. ^auffmann et al.i2003h . This bimodal population is 
composed of red, early-type galaxies with old stellar popula- 
tions and little ongoing star-formation that tend to be the most 
luminous and massive galaxies at any redshift, and a comple- 
mentary population of star-forming disk galaxies with bluer 
colors typical of young stellar populations. 

A key question is when this bimodality was established, 
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and what fraction of the total stellar mass is locked up in 
each of the two galaxy types as a function of cosmic time. 
In a landmark study. Bell et al. (2004) find that a red se- 
quence was already in place at z ~ 1. The color evolution 
of the red sequence at z < 1 is roughly consistent with pas- 
sive evolution, but analyses of the luminosity function indi- 
cate a buildup of a factor of ^2 in stellar mass over this red- 
shift range (Bell et al. 2004; Borch et al. 2006; Arnouts et al] 
lIOOTHBrown et al.ll2007t:lFaber et al . 2007; Ilbert et al.lCOlOl 
The results at z 1 are somewhat uncertain b ecause they 
requir e large corrections for incompleteness (see iFaber et aLl 
120071) . Furthermore, the most massive g alaxies do not ap- 
pear to evolve sig nificandy at < z < 1 dWake et £11120061: 
iBrown et al.l 120071) . implying that they were assembled at 
higher redshifts. 

A number of recen t studies have extend ed t his work to 
z ~ 2 In particular, lArnouts et all (12007b and lllbert et al.l 
(I20T0I) find an increase of a factor of ~ 10 in the stellar mass 
density of quiescent galaxies between z = 2 and z = 1 .2, which 
would imply very dramatic changes over a relatively short 
(w 2 Gyr) period. These studies are not definitive, as they 
depend on relatively uncertain photometric reds hifts at z > 1 
determined from broad-band NIR photometry. iTavlor et"al1 
(l2009b ) show that large redshift uncertainties make the robust 
identification of a bimodal galaxy population extremely dif- 
ficult at z > 1.5. Unfortunately, spectroscopic verification of 
these results is extremely difficult due to the faint ness of mas- 
sive g alaxies in the observer's optical (see, e.g., iKriek et akl 
l2008h . 

Here we examine the evolution of the bimodal galaxy pop- 
ulation with the NEWFIRM Medium-Band Survey (NMBS), 
a moderately deep, moderately wide near-IR survey which 
uses a novel set of medium-bandwidth filters specifically 
tuned to the redshift range 1 < z < 3 (van Dokkum et all 
120091: iWhitaker et al.1 IMO) . In iBrammer et al.l (12009b . we 
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used the NMBS to show that massive galaxies are nearly 
all red up to z ^ 2 but that a bimodal population is appar- 
ent after correcting the colors o f galaxies heavily reddened 
by dust. IWhitaker et alJ (120101) find that not only can qui- 
escent galaxies be identified in the NMBS up to z ^ 2, but 
their color scatter within the red sequence is resolved. They 
find that the fraction of quiescent galaxies among all galax- 
ies at M > 10" Mq decreases shar ply from roughly 9 0% at 
z = 1 to 40% at z = 2. Furthermore, iMarchesini et al.l (12010b 
find evidence that the quiescent fraction of extremely massive 
galaxies (M > IO^' -'M©) decreases further still to 7-30% by 
z = 3.5. 

In this paper, we use the NMBS to study the buildup of star- 
forming and quiescent galaxies from z = 2 to the present by 
constructing their mass functions and quantifying the evolu- 
tion of their number and mass densities. In ^we describe the 
survey data and sample selection. We show the galaxy rest- 
frame color distribution and its evolution in §|3] and describe 
a method of cleanly separating the red, "quiescent" galaxy 
sequence from intrinsically blue, star-forming galaxies by ac- 
counting for the effects of dust reddening. We present stellar 
mass functions in ^ We study the evolution of the galaxy 
number and mass densities in ^ and we discuss our results in 
the context of massive galaxy formation in ^ We summarize 
our results in ^ We assume a ACDM cosmology through- 
out, with n,„ = 0.3, 17a = 0.7, and Hq = 70 km s"' Mpc"'. All 
magnitudes and colors are given in the AB system. 

2. DATA 

2.1. The NEWFIRM Medium-Band Survey 

The NEWFIR M M edium-Band Survey (NMBS; 

Ivan Dokkum et all l2009t IWhitaker etaP 1201 II) provides 
well-sampled galaxy spectral energy distributions (SEDs) 
from rest-frame UV through NIR wavelengths up to z ~ 3.5, 
thanks in part to a custom set of five medium-band NIR 
filters that span observed wavelengths 1-1.7 /im at roughly 
twice the spectral resolution of standard broad-band fil- 
ters. The medium-band filter technique has been successfully 
employed at optical wavelengths to measure very precise pho- 
tometi-ic redshi fts at z 5 , 1-4 (C OMBO-17, . Wolf et al."2003|; 
COSMOS, Ilb ert et all [20091; E-CDFS, fcardamone et all 
|2010). The ability of medium-band filters to constrain pho- 
tometric redshifts depends on their sampling strong, broad 
spectral features, and in particular the Balmer/4000 A break, 
which is redshifted to A > 1 ^m at z > 1 .5. The NMBS filters 
are designed to enable precise photometric redshift estimates 
at 1.5 < z < 3.5 by improved sampling of the Balmer/4000 A 
break at these redshifts. We briefly summarize the NMBS be- 
low; a full description of the data reduction a nd photometric 
catalogs is provided bv lWhitaker et al.l (1201 II) . 

The NMBS provides medium -band near-IR pho tometry 
over ~0.25 deg^ NEWFIRM dProbst et all l200l point- 
ings in each of two well-studied survey fields, COSMOS 
(IScoviUe et al.ll2007[) and the All-wavelength Extended Groth 
Strip International Survey (AEGISO The NEWFIRM data 
were taken over 75 nights in 2008-9 at the Mayall/4m tele- 
scope at the Kitt Peak National Observatory. The NEW- 
FIRM 7123, H\2, and K data are supplemented at opti- 
mal wavelengths by deep w^r/z imaging in both fields from 
the CFHT Legacy Survejo, as reduced by the CARS team 

' ' 'http://aegis.ucoiick.org/ 
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dErben et al.ll2009l;lHildebrandt et al.ll2009l) . Furfliermore, we 
include deep Subaru imaging in hvoad-b&nd BjVjf^tz' and in 
12 medium-band filters that span 4000-8000 A (" Capak et aP 
2007; Ilbert et al. 2009). At mid-IR wavelengths, we include 
Spitzer-IKAC 3-8 /im and MIPS 24 /im imagin g that cover 
the e ntire COSMOS pointing (S-COSMOS; Sanders et dj 
120071) and 60% of the AEGIS pointing (Barmby et al. 2008). 
After masking regions with less than 30% of the maximum 
exposure time in the NMBS bands and regions around bright 
stars, the NMBS covers 0.20 deg^ and 0. 19 deg^ in COSMOS 
and AEGIS, respectively. 

Objects are detected in the NEWFIRM K image, and the 
optical/NIR photometry is performed on images convolved to 
the sam e point-spread func tion (PSF) to limit band-dependent 
effects ( Whitake r et al.l20'l 1 ; see O uadri et al.ll2007 for a sim- 
ilar photometric strategy). Objects K = 22.8 are detected at 
^5(7, and the con^esponding depths in the medium-bands are 
approximately flat as a function of f\. The CFHT and Sub- 
aru broad-band optical images are among the deepest avail- 
able in any field, which is especially important for producing 
high-S/N SEDs of optically-faint, red galaxies at z > 1 . The 
large-PSF IRAC and MIPS images require advanced photo- 
metric techniques to minimize photometric contamination by 
neighboring objects. We employ a source-fitting method that 
uses the higher-resolution /T-band image to model the posi- 
tions and sizes of objects in the redder Spitzer bands, whose 
flux normalizations are fit by least-squares regression. Fluxes 
and errors are then measured for each object with simple aper- 
ture pho tometry after subtr acti ng the model of all neighboring 
objects dLabbe et al.ll2006l; see lWuvts etani2007l for an iUus- 
trative example). The IRAC images are significantly deeper 
than the NEWFIRM /T-band. The MIPS 24 ^m images reach 
^20 ^Jy at 3(7. 

The primary sampl e used throughout this paper is defined 
as all galaxies in the IWhitaker et al.l (1201 ll) catalog (version 
5.0) of the two NMBS fields with K < 22.8, 0.4 < Zphot < 2.2, 
and the standard quality cuts, which yields 25,423 galaxies 
(15,485 with MIPS coverage). 

2.2. Photometric Redshifts 

We estimate photometric redshifts from the m-8 pm SEDs 
using the Eazy photometric redshift code (Brammer et alj 
|2008l). We use the d efault Eazy template set described by 
Bram mer et al.l (I2()08[) with a modifi ed treatment of emission 
lines inspired by lllbert et al] (120091) ; rather than using the 
emission li nes as predicted by the Pegase population syn- 
thesis code dFioc & Rocca-Volmerangdl 1997b . we compute a 
"star-formation rate" (SFR) from the rest-frame 2800 A flux 
of each template and then add Ha-j, Lyman-a, O II A3727 
and O III AA4959 ,5007 emission lines u sing fixed line ratios 
dllbert et al.ll2"009t after lKennicu"t3lI998l) . Although this treat- 
ment of emission lines is still over-simplified — real galaxies 
will have a non-trivial range of SFRs, line strengths and line 
ratios for a given 2800 A flux — we find that it significantly im- 
proves the photometric redshift quality of the medium-band 
SEDs, which are more se nsitive to line con tamination than 
broad-band SEDs (see also lllbert et al.ll2009l and Appendix lA] 
below). 

To further improve the photometric redshifts, we iterative ly 
adjust the photometric zeropoints of the ground-based photo- 
metric bands to minimize the residuals to template fits at fixed 
redshift for objects with spectroscopic redshifts (see below). 
With the exception of CFHT-m, which appears to be ~0.2 mag 
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"too faint" jErben et"ani200 9l note a similar discrepancy), the 
adjustments to the CFHT/NEWFIRM zeropoints are <2%, 
well within typical ^0.05 mag zeropoint calibration uncer- 
tainties. The Subaru broad- and medium-band images require 
significantly larger corrections to their publicly-listed zero- 
points of up to 0.2 mag (see also Table 1 of Ilbert et al. 2009). 
These offsets are constrained by the well-calibrated, overlap- 
ping CFHT photometry, and we find that the photometric red- 
shift quality is improved significantly when including all of 
the available optical data. 

Tw o large spectrosco pic surveys, DE EP2 (iDavis et al.l 
I2OO3I) and zCOSMOS (Lilly et alJ l200l . provide a large 
number of spectroscopic redshift measurements that can help 
us assess the quality of our photometric redshift estimates. 
For 2067 objects from our sam ple that have D EEP2 red- 
shifts, we measure a NMAD (see lBrammer et an i2008) scat- 
ter a/(l+z) = 0.016. In the COSMOS field, where objects 
have 35-band SEDs that include medium bands from 4000 A- 
l.T/xm, the NMAD scatter is a/(l+z) = 0.008 for 1099 galax- 
ies with zCOSMOS redshiftsQ The number of catastrophic 
redshift failures with |Az/(l+z)| > 0.1 is 2.6% (1.6%) in 
AEGIS (COSMOS). These results are among the best that 
have been achieved for photometric surveys, and we note that 
Eazy has been shown to produce redshift estimates with the 
lowest scatter and the smallest fraction of outliers of any of 
the public photometric redshift codes (Hildebrandt 2010). 

We caution that the subset of objects with spectroscopic 
redshifts tend to be optically-bright (e.g., / < 22.5 for 
zCOSMOS) and almost all are at z < 1.5; the spectroscopic 
sample is not representative of the full /iT-selected sample used 
throughout this paper van Dokkum et al. (2009) targeted four 
galaxies in the Kriek et al.l (12008) GNIRS spectroscopic sam- 
ple to test the NEWFIRM medium-band technique on K- 
selected galaxies at z > 1 .7, an d they find Az/( 1 + z) < 0.02 for 
all four objects. Furthermore. iKriek et al.l (1201 ll) have shown 
recently that Ha emission can be detected in composite SEDs 
built from the medium-band photometry itself, and they con- 
clude that the photometric redshift errors must be < 2% given 
the shape of the observed emission line. While small spectro- 
scopic samples suggest the potential of the NMBS medium- 
band filters for providing accurate photometric redshifts at 
z > 1.5, further spectroscopic followup is required to fully 
assess the redshift quality at the redshift and magnitude limits 
of the survey. 

2.3. Rest-frame colors 

We derive rest- frame U—V and y-7 color£3from the tem- 
plate that best fits the observed photometry, which itself is a 
linear combination of the Eazy templates. We take the rest- 
frame color from the template direc tly, similar to the tech- 
nique described bv lWolf et al] (l2003h for the COMBO- 1 7 sur- 
vey[3 Other methods have been developed to measure rest- 
frame fluxes by interpolating between the observed bands that 
bracket the rest-frame band at a given redshift, with the po- 
tential advantage that the resulting rest-frame fluxes are de- 
termined more directly from the observed photometry rather 
than depending on the choice of fitting templates (e.g., the 

" We use only the most reliable spectroscopic redshift quality flags from 

DEEP2 (4) and zCOSMOS (3.,v,4.,t). 

We adopt the U and V filter definitions of IMai'z ApeUanizI (2WBl . and 
2MASS-7 

EAZY2.0 with the updated templates and the rest-frame color imple- 
mentation is made available at http://www.astro.yale.edu/eazy/ 



algorithm of iRudnick et al.ll2003L imi? lemented in the INTER- 
Rest program bv lTavlor et al.ll2009bb . For a sufficiently flex- 
ible set of templates (e.g., arbitrary linear combinations of the 
Eazy templates), these two methods produce very similar re- 
sults for photometric surveys with mostly non-overlapping, 
broad-band filters. As the filters become more closely-spaced, 
direct template fluxes are more robust as they are able to make 
use of more of the observed photometry. That is, at a par- 
ticular redshift the interpolation method may interpolate be- 
tween two adjacent filters that have relatively low S/N, while 
the SED shape at the desired wavelength is actually well- 
constrained by additional nearby filters. We adopt the direct 
template U-V colors throughout, but note that our results are 
quahtatively consist ent for both our template-fitting and the 
iTavlor et al.1 (l2009bb methods. 

2.4. LiR and star formation rates 

We use the rest-frame color routine described above to de- 
rive i^I^ luminosities at 2800 A, which we combine with the 
MIPS 24 /im photometry to estimate star formation rates, 
5F7;uv+iR = 0.98 X 10"'" (Ltb + 3. 3L7snn) dBell et al.l I200I 
adapted for flie Kroupa IMF by iFranx et al.1 120081) . The 
MIPS 24/im fluxes are co nverted to to tal IR luminosities 
(LiR = Ls-iooo ^m) using the lDale & Helo u (2002) templates, 
where we adopt the log average conversion for templates with 
1 < a < 2.5 (see .Wuyts et al. 2008) . Th is conversion has a 
systematic unc e rtainty of a factor of ^3 ( Wuvts et"an i2008). 
iPapovich et"an (l2007h and lMuzzin etal J (120101) Tamong oth- 
ers) show that the commonly-u sed alternative method of fit- 
ting the luminosity-dependent IChary & ElbazI (120011) tem- 
plates tends to overestimate Ljr at z > 1 .5 by factors of ~5. 

2.5. SED modeling: stellar masses, Ay 

We measure stellar masses and dust reddening of the 
NMBS galaxies by fitting a grid of population syn thesis mod- 
els to the NMBS SEDs using the FAST code dK ri ek et"^ 
2009). The model grid is composed of iBruzual & Charlog 
(l2003i) models with exponentially-declining star formation 
histories (SF Hs) with decline rates log r/yr = 7-10, com- 
puted with a Kroupa (2001) initial mass function (IM F). To 
first order, alternative choices of the IMF (e.g. Salpeteriil955h 
cause a shift in the derived stellar masses without affect- 
ing other properties o f the fit (e.g., Msaip./MKmupa ~ 1-6; 
iMarchesini et al.l [20091) . Note that adopting alternative star- 
formation histories, such as a SFH that increases with time 
(^araston et al. 2010; Papovich et al. 2011), could also af- 
fect the derived stellar masses by significant factors, per- 
haps as large as a factor of two. We estimate the amount 
of dust reddening from the SED fit by allowing a uniform 
dust screen with up to four magnitudes of extinction in 
the V-band {Ay = 0^) and with a wavelength dependence 
following the Calze tti et alJ ([2000) reddening law. Using 
a similar SED-fitting techniq ue with photometric redshifts 
and broad-band photometry, iMuzzin et al.l (120091) demon- 
strate that st ellar masses can b e dete rmined with precision 
^0.1 dex. IMarchesini et al.l (120091) show, however, that 
changing the modeling assumptions such as the modeling 
library or the IMF can cause systematic differences in the 
derived stellar masses that are larger than these r andom er- 
rors. For example, stellar masses estimated with tMarastoiil 
( 120051) models te nd to be systemat i cally lower than those de- 
termined from Bruzual & Chariot (2003 ) models by a fac- 
tor of ~1.4 (.Wuvts et al.. .2007; .Whitaker et al.. .2010.) . We 




Figure 1. Top panels: Color-mass diagrams at 0.4 < z < 2.2. The rest-frame U -V colors are derived from the SEDs as described in 32.31 The grey lines and 
hatched regions indicate where the NMBS is less than 90% complete at the edges of the redshift bins (i.e. the region defined by the left-most of the two lines 
corresponds to the low-z edge of the bin). The number of galaxies in each bin is indicated at lower-left and the symbol size is roughly inversely proportional to 
the bin sample size for clarity. The solid red line in the top panels indicates the red-sequence found by Borch et al. (2006) at z < 1, with the redshift evolution of 
the red-sequence color zeropoint extrapolated to z = 2. Note that all massive galaxies are red, out to the highest redshifts. Bottom panels: Same as the top panels, 
but galaxies with 5FRuv+lR > 40 Mq yr"' are shown in orange. Many of the massive red galaxies at z ~ 2 are dusty star-forming systems, which means that a 
division of the galaxy population by a single rest-frame color does not lead to a homogeneous galaxy population. Note that only galaxies with MIPS coverage, 
necessary for estimating the SFR, are included in this figure. 



use the iBruzual & Charlofl (l2003h models as they appear to 
better describe the SEDs of young qui escent galaxie s than 
the Maraston (2005) models (IConrov et al. 2009: Kriek et all 
I2010I) . 

3. THE BIMODAL GALAXY POPULATION 
3.1. Color-mass relations and completeness 

We first analyze simple color-mass relations as a function 
of redshift. As we show below, the use of a single color to 
identify different galaxy populations is sufficient at low red- 
shift but is not very meaningful at redshifts z > 1 . The U -V 
color-mass relations in four redshift bins are shown in Fig. [T] 
The redshift bins are chosen such that they sample similar co- 
moving volumes, with the exception of the lowest-redshift bin 
that encompasses a volume only 40% of the other bins. 

Areas of the color-mass diagram where the sample is in- 
complete are shaded on Fig. [T] Selection in the /T-band 
is closer to a stellar mass sele ction than an optical selec- 
tion is (e.g., iFranx et aTl 120081) : nevertheless, a fixed K- 
band flux limit will result in a mass compl eteness limit 
that varies as a function of re dshift and color (iTavlor et al.l 
l2009at iMarchesini et al.l 120091) . The completeness limits 
were estimated from the NMBS sample itself co mbined with 
flie de eper FIREWORKS catalog of the CDFS (W uvts et alj 
120081) foUowing the method d escribed bv lTavlor et al. '(2009al) 
and IMarchesini et al.1 (120091) . At 7*: < 22.8 we are complete 
for all galaxies with M > lO" Mq at < z < 2.2. The com- 
pleteness limit extends to ^ 10'° Mq for blue galaxies with 



lower M/L ratios. We avoid uncertain completeness correc- 
tions below by only considering stellar mass ranges above the 
completeness threshold at a given redshift. 

The top panels of Fig. [T]highlight a key aspect of massive 
galaxies at high redshift, which has earlier been highlighted 
by others (e.g iPapovich et al.1 120061: iBrammer et all 120091 
iWhitaker et"ani2010l) : nearly aU galaxies with M > lO" M© 
have red U-V colors at 0.4 < z < 2.2. However, this does 
not imply that all these galaxies are quiescent with low star 
formation rates. In the bottom panels of Fig. [T]galaxies with 
5F/?uv+iR > 40 Mq yr"' are shown in orange. This star for- 
mation rate limit corresponds to the minimum SFR observ- 
able at z = 2.2 given the MIPS flux Hmit of 20 ^Jy. A sig- 
nificant fraction of red galaxies with M > lO^M© are vig- 
orously forming stars at z > 1, with the fraction increasing 
with redshift and reaching ^50% by z = 2 ( Whitake r et al.l 
120101) . This is qualitatively consistent with the observed in- 
crease in the number of I R-luminous galaxies with redshift 
dLe Floc'h et al.l2005Ll2009 ). In the low-z bin, galaxies on the 
massive red sequence have much lower star formation rates, 
which is consistent with many other studies that have found 
that red galaxies at z < 1 tend to be old, passi vely-evolving, 
early-type systems (e.g. ^ trateva et al. 200 il iBlanton et all 
l200llKauffmann et alJl2003l:iBell et alJI2004l) . 

3.2. Separating quiescent and star-forming galaxies 

Given that red galaxies at z > 1 comprise both quiescent 
galaxies and dusty star-forming galaxies, a single rest-frame 
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Figure 2. Dusty star forming galaxies can be cleanly separated from quiescent galaxies when using two rest-frame colors. The panels show the galaxy 
distribution in the rest-frame U — V versus V—J plane for the same redshift bins as Fig. 1. Galaxies with SFRuv+ir > 40 Mq yr"' are again shown in orange. 
These galaxies occupy a region that is dist inct from the quiesce nt galaxies. The solid line and colored bands indicate the red/blue or quiescent/dusty-Hstar-fonning 
s election developed bylLabbe et alj 120051) and lWilliams et al] t2009i) . The reddening vector for one magnitude of extinction in the V-band is indicated, assuming 
a lCalzetti et alj j2"000l) reddening law. Again, only galaxies with MIPS coverage are included in this figure. 
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Figure 3. Rest-frame U — V color distributions, with U — V corrected for the (non-evolving) slope of the color-mass relation from'Borch et al.' KOOf?) (see Fig. 
[T). The histo grams are normalized by the total number of galaxies in each redshift bin. The distribution for all galaxies with logM/MQ > 10.6 — sUghtly lower 
than our completeness limit al z = 2 — is shown in the thin, grey histograms. The top panels show the color distribution split between sources with 5Fi?uv+lR 
greater (orange, hatched) or less (black, soHd) than 4 Mp, yr~* . The bottom panels show the color distribution for red/blue galaxies selected as in Fig. [1 The 
vertical dashed line in the bottom panels indicates the lBorch et alj ^2006) red sequence selection limit, extrapolated to ^ = 2. The inset numbers show the number 
of galaxies in the "quiescent" and "star-forming" samples that have SFifuv+lR greater ("y") or less ("n") than 40 Mq yr"' . 



U-V color criterion is not the most informative way to iden- 
tify different galaxy populations. Recent studies of the lo- 
cal universe and of galaxies up to z ~ 1 have shown that the 
galaxy color bimodality is more clearly seen after account- 
ing for dust red dening (Wvderetal. 2007, Cowie & BarggB 
I2008t iMaller et al. 2009), which Brammemalj (120091) ex- 
tend to z ~ 2.2 using the NMBS. Qui escent galaxies follow 
a red sequence at lea st up to z = 2.2 dBrammer et al.ll2009t 
iWhitaker et al.l l20Tol) . while star-forming galaxies follow a 
sequence where the U-V colors become redder with increas- 
ing mass, with the color primarily determine d by increasing 
dust reddening with mass CLabbe et al..,2007l: iBrammer et al.l 

l2009h. 

iLabbe et"an (120051) show that it is possible to separate qui- 
escent from dusty galaxies with similarly red colors using 
the combination of two rest-frame colors, U-V and V-J 
(see also Wuyts et al. 2007; Williams et al. 2009). The SEDS 
of quiescent galaxies are red in U -V but blue beyond the 
Balmer/4000 A break, while dust-reddened galaxies are red in 
both colors. We choose here to use this two-color "UVJ" se- 



lection technique because it is independent of differences be- 
tween population synthesis mo dels, whi le the dust-corrected 
U-V selection used by Bra mmer et al . (2009) depends on 
accurately measuring Ay from the SED fit. jWhitaker et al.l 
(go 10) compare the UVJ and dust-corrected color selection 
methods and find that the fraction of red/quiescent galaxies 
with M > 10" Mq differs by less than 5% between them. 

Figure |2] shows the U-V and V-J colors of galaxies in 
the NMBS at 0.4 < z < 2.2. The two-color distribution is bi- 
modal up to z = 2.2 (see alsofWil liams et al. 2009; llbert et aT| 
120101) . and we adopt the lWilliams et al.l (|2009) selection crite- 
ria (solid line) to separate quiescent galaxies and star-forming 
galaxies. Galaxies with 5FiJuv+iR > 40 MQyr"' are again 
shown with orange symbols in Fig. |2] While this SFR limit, 
set by the depth of the 24 /im observations, is insufficient for 
firmly establishing a galaxy to be "quiescent", the clear sep- 
aration in the UVJ diagram supports the idea that the quies- 
cent red galaxies form a population distinct from the dusty, 
star-forming galaxies, particularly at z > 1.5. We note here 
that the existence of massive galaxies with very little ongoing 



6 



Brammer et al. 





z=[0.8,1.4] 





z=[1.4,1.8] 




z=[1. 8,2.2] 




10.5 11.0 
log M/Mg 



10.0 



10.5 11.0 
log M/Mg 



11.5 



10.0 



10.5 11.0 
log M/Mg 



11.5 



10.0 



10.5 11.0 11. 
log M/Mg 



Figure 4. Stellar mass functions for the full NMBS sample (top panels) and split using the quiescent/star-forming selection shown in Fig. [2] (middle panels). 
The points shown are simple redshift histograms divided by the volume of the NMBS, with Poisson error-bars. Representative 'Schechter' f 19761) function fits are 
shown, with the rest-frame slope fixed to a =-0.99, -1.4, and -0.7 for the full, star-fonning, and quiescent samples, respectively. The dotted lines show the local 
stellar mass function of all (black), early- (red) and late-type (blue) galaxies (B ell et al..i.20 03). scaled as described in the text. The light hatched regions show 
the 90% completeness limit for red galaxies at the high-redshift end of each bin. Note that we determine number and mass densities below by simply counting 
objects at masses where the NMBS is complete, rather than integrating the Schechter functions. The bottom panels show the fraction of red, quiescent galaxies as 
a function of stellar mass and redshift. The dotted line shows the ratio of the Bell et al. 1 2003) early and early-l-late Schechter functions, while the solid lines show 
the ratio of the Schechter function fits to the NMBS mass functions. Quiescent galaxies cleai'ly evolve much more rapidly than star-fonning galaxies, driving the 
evolution of the total mass function at the high mass end. Here and in the figures below, all galaxies are shown whether or not they have MIPS coverage. 



star-fomiation at z 2 is supported by near-IR spectroscopic 
studies (Kriek et al. 2006). 

It is clear that red (quiescent) and blue (star-forming) sam- 
ples selected using the UVJ colors will be different from 
those based on a sing le rest-frame color, as was done by, e.g., 
iBellet all ( 12004 and lBorch eta[\ (l2006h . This is shown qual- 
itatively in Fig.|3] which shows histograms ofU-V color cor- 
rected for the slope of the color-magnitude relation. The top 
panels show the color distribution split according to the SFR 
threshold used in Fig. [T] while the bottom panels show the 
distribution split according to the UVJ quiescent/star-forming 
selection. Although the distribution of U -V colors changes 
only slightly with redshift, the distribution of star-forming 
versus quiescent galaxies changes rapidly with redshift. Here- 
after, we discuss the two distinct populations apparent in the 
UVJ diagram as "quiescent" and "star-forming" samples. We 
cannot exclude some on-going star formation in the quiescent 
galaxies (particularly if it is completely obscured), but their 
rest-frame optical SEDs are dominated by an evolved stellar 
population. 

It is worth commenting here on other recent studies of the 
evolution of the red sequence to z ~ 2. A number of stud- 
ies have found little or no evidence fo r a significant num- 
ber of red, quiescent galaxies at z > 1.5 (lArnouts et al.ll200"7t 
ICirasuolo et al."2007). Like the NMBS, these studies rely on 
photometric redshifts, but they have relatively poor sampling 
of observed NIR wavelengths with only the / and K broad- 
band filter s. This will result in l arge redshift uncertainties at 
z = 1.5-2 dBrammer et al.l 120081) . and therefore derived rest- 
frame colors will have insufficient precision to ident ify the 
distinct red/blue populations at those redshifts (Taylor et alJ 



l2009ah . Though sample sizes are small, spectroscopic stud- 
ies have unambiguously identified quiescent galaxies at z 2 
(ICassata et a l. 2008; Kriek et al. 2008), emphasizing the need 
for precise redshift measurements. Furthermore, the def- 
inition of "quiescen ce" varies from one study to another 
lArnouts et al.l (120071) define galaxies to be quiescent when 
they are best-fit by a non-evolving (local) elliptical galaxy 
template. As can be seen in Fig.[T] passive evolution ensures 
that quiescent galaxies at z > 1.5 cannot have the colors of 
local elliptical galaxies, and thus such a selection will fail to 
identify this population. 

4. STELLAR MASS FUNCTIONS 

Having robustly divided the sample in quiescent and star- 
forming galaxies, we now study the evolution of the stel- 
lar mass function split by galaxy type. We first show 
the evolution of all galaxies, in the top panels of Fig. [4] 
(see also, e.g., Font ana et alj l2006t iMarchesini et all 12009^ 
Pozzetti et al. 2010). These mass functions are computed by 
simply counting galaxies in stellar mass and redshift bins. We 
do not adopt the Vmax formalism (.Avni & Bahcall 1980) as we 
only c onsider stellar m asses where the NMBS is complete. 
We fit ISchechte^ (1197 6*) functions with fixed faint-end slope 
to the densities in each redshift bin to demonstrate only that 
the mass functions have reasonable shapes. Our stellar mass 
completeness limits do not allow us to constrain the faint-end 
slope at z > 1, a nd strong degenera cies between the Schechter 
parameters (e.g. IMarchesini et al.ll2b09.) would make the pa- 
rameter values and their evolution with redshift difficult to in- 
terpret given the simple analysis used here. We defer a more 
detailed analysis of the NMBS stellar mass functions, includ- 



Galaxy number and mass densities to z = 2.2 



7 



10' 



10" 



Q. 



10' 



Apr 1 3 1 0:1 1 ;07 201 1 ] 
I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' 

□ o 

□ 

Y7 V ° ■& ® ^ 1 



m Taylor et al. 2009 
V PozzeftI et al. 2010 



O 




0.0 



0.5 



1.0 



1.5 



2.0 



2.5 



8.5 - 



8.0 - 



7.5 - 



o 



7.0 - 



i □ 



i 



X ^ ^ 



■ u — 



X Boroh et al. 2006 
□ llbert et al. 2010 
<> Maroheslnl et al. 2009 
« Bell et al. 2003 



o 




_1_ 



_1_ 



0.0 



0.5 



1.0 



1.5 



2.0 



2.5 



Figure 5. (Left panel) Number density evolution of galaxies in tlie NMBS selected to mass limits of logM/Mg > 10,11 (large filled circles with eiTor bars). 
The redshift error-bars show the width of the redshift bin used, which are half the width of those used in Figs[l]2] The density eiTors include Poisson and cosmic 
variance errors computed following Somerville et al. (2004), added in quadrature. The densities for the individual AEGIS and COSMOS fields are shown with 
dotted and dash-dotted lines, respectively. (Right panel) Mass density integrated to the mass limits indicated in the left panel. The eiTors shown are the fractional 
errors on the number densities, and do not include systematic uncertainties associated with the stellar-mass determinations. Number and mass densities integrated 
from the mass functions of Borch et al. (2006); Marchesini et al. (2009); llbert et al. ( 2010), scaled to match our SED modeling assumptions where necessary, are 
shown with symbols as labeled. Number densities are taken directly from Tavlor et al., ( ;2009a) and .Pozzetti et al.. ( .2010.) . The numerical values of the densities 
and total errors for the full NMBS sample are provided in Table[T]below. 



ing a full accounting of systematic errors and incompleteness, 
to a future paper 

We inc l ude in Fig. |4] the z = 0.1 mass function from 
iBell et all (|2003|) for comparison, which we have scaled to 
our as sumed cosmology. Additionally, we scale the IBell et al.l 
(|2003') stellar masses down by a factor of 1.2 to account for 
the difference between a "diet-Salpeter" IMF and the Kroug^ 
(|2001) IMF we use to estimate stellar masses. It is apparent 
from the total mass functions that the mass function evolves 
gradually from z = to z = 2, with no indications for sudden 
or dramatic changes in particular redshift ranges. 

The mass functions are split into quiescent and star-forming 
galaxies in the middle panel s of Fig. |4] The dotted lines 
indicate the iBell et all (l2003h local mass functions of early- 
and late-type populations separated by color It is clear that 
the mass functions of the quiescent and star-forming sam- 
ples evolve in significantly different ways. Massive galaxies 
(M > 1O"M0) at z < 1 are almost entirely in the quiescent 
population. The number of massive quiescent galaxies de- 
creases steadily with increasing redshift, while the mass func- 
tion of the star-forming sample evolves very little up to z = 2 
forM > 10""^ Mq. The differential evolution is such that the 
two populations are approximately equal in number at z = 2 
for these galaxies with M > 10'° "^Mq (Fig. |4] bottom pan- 
els). Thus we obs erve in the galax y population directly the 
argument made by ' Bell et all (120071) . who noted that the IR- 
luminous (massive) starburst galaxies observed at z > 2 must 
later migrate to the red sequence to avoid dramatically over- 
predicting the number of massive, blue galaxies observed lo- 
cally. Furthermore, we demonstrate here that the evolution of 
the total galaxy population is mostly due to the evolution of 
the quiescent galaxies: the star-forming population is similar 
at all redshifts, and it is the rise of quiescent population with 
cosmic time that is responsible for the gradual evolution of 
the total mass function. 

The shape of the quiescent mass function may also evolve 
with redshift, although we are severely hampered by incom- 



pleteness at low masses and high redshifts. Taking the data 
at face value, we find that massive quiescent galaxies above 
10"Mo evolve more slowly than those at lower masses, but 
deeper data at the low mass end and spectroscopic redshifts 
at the high mass end are needed to confirm this. We note 
that this evolution would be qualitatively consistent with the 
recent (z < 1) build-up of low-mass red galaxies found by 
iRudnick et al.l (l2009h . both in cluster and field populations. 

5. NUMBER AND MASS DENSITY EVOLUTION 
5.1. Full galaxy sample 

To further quantify the evolution of star-forming and quies- 
cent galaxies, we now consider the evolution of their number 
and mass densities with redshift. We measure these quantities 
directly from the data by simply counting objects above the 
mass completeness threshold in a given redshift bin, so we do 
not rely on the assumption that the mass functions follow a 
Schechter function. Again, the NMBS does not constrain the 
faint-end of the mass function at z > 1, so integrating the mass 
functions to zero mass would result in an unreliable extrapo- 
lation of the observed data. 

The number and mass density evolution in two mass ranges 
is shown in Fig.|5] The redshift bins used are those from Figs. 
[T]-0] divided in two; the bins are still somewhat larger than 
t he ex pected photometric redshift uncertainties even at z = 2 
O2.2I 1. Only redshift bins where the survey is complete to 
the specified mass limits are shown. The uncertainties on the 
number densities are a combination of Poisson errors and the 
cosmic variance estimated using the ISomerville et al.l (120041) 
prescription, given the observed number density, redshift bin- 
ning, and survey geometry. We can obtain a rough estimate 
of the cosmic variance directly from the NMBS itself: the 
number and mass densities for the individual COSMOS and 
AEGIS fields are shown as the dotted and dash-dotted lines, 
respectively. The differences between the two fields are sig- 
nificant and are generally well-represented by the error esti- 
mates that account for cosmic variance. 
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Figure 6. Number and mass density evolution for galaxies divided into quiescent and star-forming samples as in Fig. |2]for three stellar mass bins. Redshi ft bins 
are only shown where the mass bins ai'e > 90% complete. The z = 0.l density measurements ai'e integrated from the early- and late-tvpe lBell et afl 120031) mass 
functions (see Fig. [3] Number densities of the "red" and "blue" types defined by Pozzetti et al. 1 2010) are shown with open triangles, taken from their Fig. 13. 



Note that the Pozzetti et al.i 12010 ) mass bins are slightly different than those used here (lO'" 
and total errors for the full NMBS sample, divided by galaxy type, ai'e provided in Table|2]below. 
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The inclusion of the local density measurement is impor- 
tant for evaluating the overa ll evolution , and here we include 
densities integrated from the lBell et al.l (f2003) mass function. 
While the methodology we use to estimate stellar masses and 
separate quiescent fro m star-form ing galaxies is quite differ- 
ent than that used by iBell et all (|2003), the local reference 
point appears to connect well with our NMBS measurements 
in Fig. |5] Similar to other studies, we find that there is 
modest mass density evolution (0.25 dex) up to z = 1 for 
masses M > 1O'°M0. Up to z = 2, we observe an overall 
decrease in the mass density of ^0.6 dex at M > lO" Mq. 
We do not find a sudden change at z > 1; the evolution ap- 
pears to be gradual over the entire redshift range < z < 2.2. 
We note t hat the resul ts here are consistent with those of 
iMarchesi ni et al.l (120091) . who find very little evolution at ex- 
treme masses > lO"'^ Mq (Figure|4|l. 

Differences between the two NMBS fields highlight the 
importance of sampling large volumes: the COSMOS field 
shows an overall mass density decrease of 0.6 dex between 
0.4 < z < 2, while densities measured in the AEGIS field are 
essentially consistent with no evolution over this same red- 
shift range. Furthermore, we note that our results are sensitive 
to subtle redshift-dependent systematic errors in the masses, 
including systematic differences between the z = point and 
the higher redshift data. This is demonstrated explicitly in 
Fig.[TO]in the Appendix. 

5.2. The number and mass densities of quiescent and 
star-forming galaxies 

We examine the number and mass density evolution of the 
quiescent and star-forming samples ( ^3.2t in Fig. |6] Rather 
than integrating down to a mass limit as in Fig. |5] we now 



consider three separate stellar mass bins where the NMBS is 
complete to at least z = 1 . While this is essentially the same 
measurement as the stellar mass functions described above, 
plotting the densities allows the trends with redshift and mass 
to be more readily apparent. 

It is immediately clear again that the quiescent and star- 
forming galaxy populations evolve in very different ways. 
The number and mass densities of star-forming galaxies re- 
main nearly constant with redshift for all masses log M/M0 > 
10.2. Star-forming galaxies with logM/M0 > 11, which are 
relatively rare at low redshift, are more prominent at z = 2 
where they have the same number and mass densities as the 
massive qui escent population . We again take the local com- 
parison from lBell et al.l (120031) . who provide stellar mass func- 
tions divided into early- and late-type populations by color. 
The local comparison is important for showing the overall 
redshift evolution, particularly at low masses where our data 
indicate rapid growth of the quiescent population since z ^ 1 . 
Whether or not the local comparison is considered, the density 
evolution of star-forming galaxies is in stark contrast with that 
of quiescent galaxies, whose density decrease with redshift at 
all of the masses considered. The evolution of the mass den- 
sity of the quiescent galaxies as a function of redshift is to first 
order independent of mass for \ogM /Mq > 10.2, with about 
0.5 dex of mass density growth per unit redshift for all of the 
mass ranges shown. Because three times as much time elapses 
between z = 1 and the present day as between z = 2 and z = 1, 
the massive red sequence grows rapidly between z = 2 and 
z = 1 , at which point the growth slows considera bly over the 
remaining ~8 Gyr until the present day (see also iKriek et aLl 
I2OO8I) . 

The differential density evolution between the quiescent 



Galaxy number and mass densities to z = 2.2 



9 



T 



log n(M>M.) [Mpo"'] 




I I 1 1 I I I I I I J 

1 2 

z 

Figure 7. Redshift evolution of stellar mass at fixed cumulative number den- 
sity, n{M > Mc), for the quiescent and star-forming galaxy samples. The un- 
certainty on Mc, estimated from Monte-Carlo simulations in which we per- 
turb the stellar masses with a random error of 0.1 dex and remeasure the 
derived quantity, is indicated by the width of each col or/density track. The 
tracks are extended to z = using the lBell et alj <20Q3I) mass functions. The 
shaded region at z=0 demonstrates the extent of the systematic uncertainty 
related to the model library used to derive stellar masses. The gray shaded 
region shows the completeness limit of the NMBS. 

and star-forming populations shown in Fig. |6] provides an- 
other manifestation of the "downsizing" phenomenon, in 
which the characteristic sites of star-forrnation s hift to lower 
mass galaxies at later times ^Cowie et al."'1996", followed by 
e.g., Ifelletal. 2005; Juneau et al. 2005; Bundv et al.l 120061 
and many others). Here we see that the number and mass den- 
sities of quiescent and star-forming galaxies are equal at ear- 
lier times for increasing galaxy stellar mass, after which the 
densities are dominated by the quiescent population. This is 
effectively a reformulation of taking the stellar mass at which 
the mass functions cross to be the characteristic mass of trans- 
forming gal axies from t he star-forming to quiescent popula- 
tion (Bund y et alJl2006h . Indeed, we see a similar trend in 
the NMBS stellar mass functions in Fig. HI the quiescent and 
star-forming stellar mass functions cross at stellar masses that 
increase with redshift. 

6. DISCUSSION: THE BUILD-UP OF QUIESCENT GALAXIES 

The trends for star-forming and quiescent galaxies are con- 
trary to naive expectations. One would naively expect that 
the stellar mass density of star-forming galaxies increases 
with time, as they are forming new stars. Quiescent galax- 
ies have already stopped forming stars, and one might expect 
their mass density to remain constant. However, we observe 
the o pposite. Consistent with previous s tudies at lower red - 
shift ( Borchet al.ll2006t iBellet al. 1 120071: iMartin et al.ll2007h . 
we find that the mass density of the star-forming popula- 
tion shows very little evolution up to z 2 for galaxies with 
logM/M0 ^ 10.6. At the same time, the mass density in qui- 
escent galaxies increases with time. The obvious interpreta- 
tion is that galaxies migrate from the star-fo rming population 
to the quiescent population (see also, e.g., iBell et all 120071 ; 
iFaber et aDl2007 ). Furthermore, the mass density of the mas- 
sive quiescent population above 1O"M0 might grow through 



mergers with less massive galaxies (e.g.. Ivan Dokkuml2005l) . 

6.1. Mass evolution affixed cumulative number density 

To explore the causes of the rise of quiescent galaxies we 
first ask how much the masses of individual gala xies grow 
with time. As shown by Ivan Dokkum et al] (120 IQ) . select- 
ing galaxies at a constant number density rather than a con- 
stant mass enables the study of the evolution with time of a 
single coherent population of galaxies at the massive, expo- 
nential end of the mass function, as both star formation and 
any merger with a ratio somewhat less than than 1 : 1 will in- 
crease th ese galaxies' masses witho ut changing their number 
density. Ivan Dokkum et al] (120101) find that galaxies with a 
number density « = 2 x 10""* Mpc"^ dex"' grow in mass by a 
factor of two sinc e z = 2 (with logM/ Mii^ = 11.15 atz = 2). 

We modify the Ivan Dokkum et al.l (1201 Oh selection slightly 
and consider a selection based on a constant cumulative num- 
ber density, defining Mc to be the stellar mass at constant 
n{M > Mc). This definition has the advantage that it is single- 
valued even at low-masses. Starting with a baseline mass 
function with Mc = Mc,o at time t = 0, the value of Mc can 
increase after time t (i.e., M^,, > Mc.o) as a result of three phys- 
ical process that all may occur simultaneously: growth from 
mergers where the secondary objects have M < Mc or both 
progenitors have M w Mc, star formation, and creation of new 
galaxies in the considered sample via transformation of galax- 
ies (with M > Mc) from a different population. In the first two 
of these processes, an increase in Mc,, is directly related to 
the average increase in mass in individual galaxies above that 
threshold, given the simplifying assumption that all galaxies 
above the threshold evolve in a uniform way. If two galaxies 
each with M > McX) merge within time, t, then Mc,, < Mc.o, 
as one would have to move further down the mass function to 
count the same number of galaxies in the same volume. Such 
mergers with both progenitors with M > McX) would necessar- 
ily have mass ratios of nearly 1 : 1 for Mc 10" Mq due to the 
steepness of the mass function. Finally, mass loss in galaxies 
above McX) after some time, t, will result in Mc,, < Mc.o- 

We show the evolution of Mc for the quiescent and star- 
forming galaxy populations in Fig. |7] for a range of n{M > 
Mc). We find that the stellar mass of quiescent galaxies with 
log M/M0 > 1 1 at z = 2 in creases by a factor of ^2 by z = 0, 
similar to the lvan Dokkum et a l. (2010) value for the full sam- 
ple not divided by color. Interestingly, the growth of Mc at 
the highest masses slows with time, perhaps suggesting an in- 
creased prevalence of dry, equal-mass mergers at these masses 
since z = 1. The value of Mc for the star-forming galaxies 
is nearly constant with redshift for all densities considered, 
which is likely the result of tension between competing and 
opposite effects: star-formation within the population and mi- 
gration out of the star-forming and into the quiescent popula- 
tion. 

6.2. Mergers versus truncation of star formation 

We now construct a very simple empirical calculation to test 
the hypothesis that the observed number density evolution of 
quiescent galaxies is a result of growing the stellar mass of 
individual galaxies by the factor of two measured from Fig. |7] 
For the quiescent population, this growth should be dominated 
by minor mergers since 1) these galaxies are not forming stars 
at a high rate and 2) equal mass mergers are rare, at least 
at the massive end, because the galaxies themselves are rel- 
atively rare (see also .van Dokkum et al...2010.) . The shape of 
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Figure 8. Number density evolution of quiescent galaxies. The shaded regi ons are the same number densities and their uncertainties for quiescent galaxies as 
in Fig. |6] The stars indicate the z = 0. 1 measurement from ' Bell et al.l 00031) . The solid lines indicate the num ber d ensity evolution produced by a very simple 
calculation in which we assume that all quiescent galaxies grow in mass by a factor of two from z = 2 to z = ( j|6.1) . 

the mass function ensures that there i s a large number of low- 
mass merger progenitors. Taking the iMarchesini et all (120091) 
1 .3 < z < 2 function, we measure the change in number den- 
sity in each mass bin that results from increasing all masses 
by a factor of two. This is essentially taking the derivative of 
the mass function at the indicated mass bins. Thus, the num- 
ber density evolution for a given increase in stellar mass is 
greatest in the mass bin where the mass function is steepest. 

The density evolution that is implied by this calculation 
is shown in Fig. [8] where we normalized the evolution the 
density observed in the NMBS at z = 1. The modest as- 
sumed growth of stellar mass in individual galaxies results 
in a factor of 3 growth in the number density since z = 2 
for the massive bin. This simple model can explain at least 
^50% of the observed density growth in the massive bin, 
suggesting that the growth in the number and mass densities 
of these ga laxies is d ominated by (minor) m ergers (see also 
iHopkins et al. 2010; van Dokkum et alj2010l) . The remaining 
density growth of massive quiescent galaxies can likely be ex- 
plained by transforming galaxies from the dusty star-forming 
population, which is indeed required to avoid overproducing 
massive star-forming galaxies, as in Fig. |4] 

By contrast, growing the stellar masses of individual quies- 
cent galaxies by a factor of two can only account for 10-15% 
of their observed number density evolution at lower masses. 
If this modest mass growth within the less massive quiescent 
population (e.g. due to red mergers) is insufficient to explain 
the rapid increase of their number and mass densities, trans- 
formation of galaxies from the star-forming to quiescent pop- 
ulations is likely the dominant mechanism necessary to make 
up the difference. This conclusion is further supported by the 
lack of evolution of the number/mass densities and for the 
star-forming population as seen in Figs. |6]and|7] respectively. 

6.3. Consequences of the massive, dusty galaxy population at 

z>l 

The process by which galaxies transform from the star- 
forming sequence to the quiescent sequence is unknown, but 
at moderate and high stellar masses it can not be a simple 
fading of blue, dust-free star forming galaxies to red, dust- 
free quiescent galaxies. As noted earlier, the bimodal galaxy 
population at z > 1 is not immediately apparent in the U-V 



color-mass diagram because the star-forming and quiescent 
galaxies have essentially identical red colors (Fig. [T]i. The 
former, we have argued, are heavily reddened by dust. At no 
point for z < 2 do we see a significant population of massive, 
star-forming galaxies with un-reddened blue colors. This lack 
of observed blue galaxies indicates that the dusty nature of 
the massive star-forming galaxies we observe at z > 1 is not 
likely to be a simple effect of the dust geometry — for example 
that highly-inclined (blue) disks would have reddened colors 
that would put thern_onJ;he red sequence (e.g. Bailin & Harrii 
2008; Mail er et al.i2(j09l) — ^because the random orientation of 
objects would show a range of obscuration at constant stellar 
mass broader than that observed. 

As a result, massive galaxies must leave the star-forming 
sequence and eventually arrive on the quiescent sequence at 
essentially the same U-V color Although local early-type 
galaxie s are frequently observed to co ntain modest amounts 
of dust dKormendv & Diorgovskill989[) . the dust does not ap- 
pear to affect the overall optical colors of nearby red-sequence 
galaxies to the dramatic extent that we find in the optical- 
NIR SEDs of massive star-forming galaxies at z > 1 (see also 
Whitakeret al. 2010). Therefore, dust must be destroyed or 
the dust geometry needs to change during the process of of 
transforming galaxies from the dusty star-forming sequence 
to the quiescent red sequence. As these processes probably 
operate on timescales similar to or longer than the timescale 
of star formation quenching, one might expect an increase in 
the amount of du st in the quiescent population with redshift. 
IMarchesini et al. (2010) do find that extremely massive galax- 
ies at z > 3 in the NMBS have fairly strong Balmer/4000 A 
breaks and also > 1 mag of dust extinction, measured with 
the same techniques described in 32.51 Future work combin- 
ing medium-band surveys such as the NMBS with new far-lR 
capa bilities to better ch aracterize emission from the dust itself 
(e.g. lElbaz et aUlTolOl) . will help to further test this prediction 
of moderate dust in quiescent galaxies at z > 1 . 

7. SUMMARY 

We use the unique photometric dataset provided by the 
NEWFIRM Medium-Band Survey to study the color and stel- 
lar mass distributions of galaxies up to z = 2. The NMBS 
provides the best-sampled SEDs and most precise photomet- 
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ric redshift estimates at z > 1 available to date, which greatly 
improves the quality of rest-frame properties (colors, masses, 
stellar population parameters) derived from the photometry. 
We find that nearly all galaxies with logM/M© > 10.5 have 
red rest-frame U -V colors up to z = 2. Using a two-color 
method, supported by MIPS 24/im photometry, we identify a 
bimodal galaxy population consisting of a quiescent sequence 
with relatively low star-formation rates and a distinct star- 
forming sequence that becomes increasingly dusty with in- 
creasing stellar mass. 

Separating the evolution of the mass function by star for- 
mation rate, we find that it is driven by the rise of quiescent 
galaxies from z = 2 to z = 0. The mass function of star-forming 
galaxies is remarkably similar at all redshifts, whereas the qui- 
escent galaxies show strong, mass-dependent, evolution (see 
also [Marchesini et al. 2009; Droryetal. 2009). Quantifying 
this evolution, we find that number and stellar mass densi- 
ties of all galaxies with logM/M0 > 10 evolve by ^-^0.5 dex 
per unit redshift. Considering the density evolution of star- 
forming and quiescent galaxies separately, we find that the 
density of star-forming galaxies is nearly flat out to z = 2, and 
the density of quiescent galaxies decreases by a factor of 10 
from z = to z = 2. Interestingly, at z = 2 the mass and number 
densities of the quiescent and star-forming populations with 
M > 10" Mq are nearly identical. This crossing point of the 
densities occurs at progressively later times at lower stellar 
masses. 

Using an empirical argument based on selecting subsam- 
ples based on their cumulative number density, we show that 
the average mass in individual quiescent galaxies grows by 
a factor of ~2 from z = 2 to z = 0. If we assume that most 
of this mass growth is due to mergers, we find that a sim- 
ple calculation based on the shape of the mass function is 
able to explain much (at least ^50%) of the number den- 
sity evolution of galaxies with M > 1O"M0. This model is 
unable to account for the density evolution of less massive 
quiescent galaxies, which we argue are formed primarily via 
transformations from the star-forming population. In general, 
our results bring together a large variety of results from the 
literature — ^based on disparate samples selected at different 
wavelengths and redshifts — into a coherent picture extended 
up to z = 2, thus demonstrating the utility of large, uniformly- 
selected surveys for the study of galaxy formation and evolu- 
tion. 
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Figure 9. Emission line contamination of broad/medium-band fluxes. The broad-band photometry of each object with 0.4 < z < 2.2 and logM/MQ > 9.5 
are shifted into the rest-frame and the SEDs are normalized by subtracting a linear fit to the computed rest-frame V and i'-band fluxes. Star-forming (blue) and 
quiescent (red) galaxies are selected as in Fig. |2] and the shaded areas show the l-cr range of the observed fluxes shifted into the rest-frame. An Ha emission 
feature is clearly visible for the galaxies selected to be star-fonning, even though this feature is only crudely sampled by the i-K broad- and medium-band filters 
over the redshift range shown. The thin red and blue lines show two Eazy templates integrated through the NEWFIRM J2 filter at < z < 4. 
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ARPENDIX 

A. LINE EMISSION IN THE BROAD/MEDIUM-BAND FILTERS 

We find that adjusting the treatment of nebular emission lines in the Eazy templates, following lllbert et aTl (12010 1). significantly 
improves the quality of our photometric redshift estimates. This treatment becomes even more important as we use medium-band 
photometry from 4000A to 1.7/im and the medium-width filters are more sensitive to the flux contamination by emission lines. 
We investigate the contribution of emission lines to the medium- and broad-band photometry in Fig. |9] The shaded regions show 
the distributions of b road-band fluxe s shifted to the rest-frame of objects in the quiescent and star-forming populations, selected 
as in FigfTICsee Whitaker et al.l ( 120101) for a more detailed description of the construction of these "rest-frame SEDs".) We include 
galaxies over the full redshift range of the NMBS and logM/M© > 9.5. We detect the clear signature of Ha emission, with a 
shape similar to that expected from a typical filter sensitivity curve, in the average SED of the star-forming sample. This feature 
is sampled in the figure by, among others, the Ji, J3, and H2 filters at z = 0.6,0.9, and 1.6, respectively. The quiescent population 
shows no such feature and has a different overall SED shape from the star-forming population. These results highlight both 
the importance of including emission lines in the redshift templates and also that our star-forming/quiescent selection criterion 
succeeds in selecting galaxies with/without Ha emission. The detection of spectral features in the NMBS photometr y, including 
the shape and strength of the Ha emission line as a function of galaxy SED type, is explored in greater detail by iKriek et alj 
(12011.) . 

B. EFFECT OF MASS MEASUREMENT ERRORS ON THE DENSITY EVOLUTION 

Redshift-dependent systematic errors in the stellar masses are an important systematic uncertainty in our study. It is very 
difficult to assess how important such errors are. As noted in, e.g. .'March esini et al.l (120091) the metallicities, IMF, and other pa- 
rameters could vary systematically with redshift. Furthermore, we sample a different part of "model space" at low redshift (where 
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Figure 10. Effect of mass measurement en'ors on number and mass density estimates. Tlie shaded areas indicate the range of densities measured adopting the 
indicated mass cuts plus or minus a systematic measurement "error" of 0. 1 dex. 



Table 1 

Number and mass densities 



IogM/M0 > 10 IogM/M0 > 11 

z n" n p 



0.61 34.4±5.2 13.9±2.1 2.5±0.6 3.8±0.9 

0.96 31.0±4.7 13.2±2.0 2.2±0.6 3.5±0.9 

1.23 ■■■ ■■■ 1.4±0.3 2.1±0.5 

1.48 ■■■ ■■■ 1.7±0.4 2.7±0.6 

1.69 ■■■ ■■■ 1.6±0.4 2.6±0.6 

1.90 ■•■ ■■■ 1.5±0.3 2.4±0.6 

2.10 ■■■ ■■■ 1.2±0.3 2.0±0.5 

2.28 ■■■ ■■■ 0.7±0.2 1.2±0.3 



" n : Mpc-^ p: lO' Mq Mpc^^ 



Table 2 

Number and mass densities, split by galaxy type 



10.2 <logM/M0 < 10.6 
rf p" 



10.6 < logM/Mg < 11.0 
n p 



11.0 < logM/Mg < 11.6 
n p 



z 


Q" 


SF'' 


Q 


SF 


Q 


SF 


Q 


SF 


Q 


SF 


Q 


SF 


0.613 


7.2± 1.4 


7.3± 1.4 


1.9±0.4 


1.8± 0.3 


6.5± 1.3 


2.1±0.5 


4.1±0.8 


1.2± 0.3 


2.2± 0.5 


0.3±0.1 


3.4± 0.8 


0.3± 0.2 


0.961 


6.0± 1.2 


6.3± 1.3 


1.6± 0.3 


1.6± 0.3 


5.9± 1.2 


3.1±0.7 


3.8±0.8 


1.8± 0.4 


1.9±0.5 


0.3±0.1 


3.0± 0.8 


0.4± 0.2 


1.234 


3.1± 0.6 


6.3± 1.0 


0.8±0.2 


1.6± 0.3 


3.1±0.6 


2.1±0.4 


2.0± 0.4 


1.3± 0.3 


1.0± 0.3 


0.4± 0.1 


1.5±0.4 


0.5±0.2 


1.477 










2.9±0.5 


2.4± 0.5 


1.8±0.3 


1.4± 0.3 


1.2±0.3 


0.5±0.1 


2.0± 0.5 


0.7± 0.2 


1.692 










2.4± 0.5 


2.7±0.5 


1.5±0.3 


1.6± 0.3 


1.0± 0.3 


0.6± 0.2 


1.7±0.4 


0.8±0.3 


1.900 


















0.7± 0.2 


0.8±0.2 


1.1±0.3 


1.2±0.3 


2.100 


















0.6± 0.2 


0.6± 0.2 


0.9± 0.3 


0.9±0.3 


2.275 


















0.2± 0.1 


0.4± 0.1 


0.3±0.1 


0.7± 0.2 


2.425 


















0.2±0.1 


0.5± 0.2 


0.3±0.1 


0.9± 0.3 



" n: \Q-* Mpc-3, p- 10^ Mq Mpc"' 
Q: quiescent, SF: star-fonning (Fig. |2). 

we rely on models for 10+ Gyr old stellar populations) and high redshift (where we rely on models for younger populations, and 
star-forming galaxies). In Fig. [TO]we show the effects of a 0.1 dex error in the masses on the number densities and the mass 
densities. The effects are substantial, particularly for the number densities at the highest masses. A 0.1 dex redshift-dependent 
error implies a range in the evolution of the number density of massive galaxies from very little to a decrease of a factor of ^5 
between z = 0.4 and z = 2. The only way to reduce these uncertainties is to obtain dynamical mass measurements of large samples 
of high redshift galaxies. 



